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ABSTRACT 

We present detailed optical, near-infrared, and radio observations of the X-ray 
flash 050416a obtained with Palomar and Siding Springs Observatories as well as 
the Hubble Space Telescope and Very Large Array, placing this event among the 
best-studied X-ray flashes to date. In addition, we present an optical spectrum 
from Keck LRIS from which we measure the redshift of the burst, z = 0.6528. At 
this redshift the isotropic-equivalent prompt energy release was about 10 51 erg, 
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and using a standard afterglow synchrotron model we find that the blastwave 
kinetic energy is a factor of 10 larger, E^^so ~ 10 52 erg. The lack of an observed 
jet break to t ~ 20 days indicates that the opening angle is 6j > 7° and the total 
beaming-corrected relativistic energy is > 10 50 erg. We further show that the 
burst produced a strong radio flare at t ~ 40 days accompanied by an observed 
flattening in the X-ray band which we attribute to an abrupt circumburst density 
jump or an episode of energy injection (either from a refreshed shock or off-axis 
ejecta). Late-time observations with the Hubble Space Telescope show evidence 
for an associated supernova with peak optical luminosity roughly comparable 
to that of SN1998bw. Next, we show that the host galaxy of XRF 050416a 
is actively forming stars at a rate of at least 2 M yr _1 with a luminosity of 
Lb ~ 0.5L* and metallicity of Z ~ 0.2 — 0.8 Z & . Finally, we discuss the nature 
of XRF 050416a in the context of short-hard gamma-ray bursts and under the 
framework of off-axis and dirty fireball models for X-ray flashes. 

Subject headings: gamma-ray bursts: specific (XRF 050416a) 



Introduction 



Nearly a decade ago, X-ray flashes (XRFs) were observationally recognized as a subclass 
within th e sample of gamm a-ray bursts (GRBs) detected by the BeppoSAX Wide-Field 
Cameras (IHeise et al. Il200ll ). The events are distinguished by a prompt spectrum that peaks 
in the soft X-ray range (E p < 25 keV), a factor of ~ 10 below the typical values observed 



for GRBs (IBand et al. Ill993l ). Since then, it has been shown that XRFs and GRBs shar e 



many obs ervational properties, inc luding prompt emission durations (ISakamoto et al. 



2005 



redshifts (ISoderberg et al. Il2004af) . broadband afterglo ws (e.g., XRF 050406. iRomano et al 



20061; 


Schadv et al. 


2006), and host 


2004 


Rau. Salvato & Greiner 


2005) 



(SNe) in association with XRFs 020903 JSoderberg et al. I I2005M iBersier et al. 1 120061) and 



060218 flPian et al. 1120061 : Mirabal et al. 1120061 : Modiaz et al. 1120061 ) indicate that XRFs, like 



GRBs, are produced in massive stellar explosions. Together, these clues strongly suggest that 
XRFs and GRBs share similar progenitors. 

Driven by this progress, several theories have been proposed to explain the soft prompt 
emission observed for XRFs under the framework of a standard GRB model. One popular 
idea posits that XRFs are merely typica l GRBs viewed away from the collimation axis (e.g., 
Yamazaki. Yonetoku fc Nakamural 120031 ) . In this scenario the prompt emission is primarily 
beamed away from our line-of-sight, resulting in lower fluence and E p values for the observed 
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burst. An important implication of the off- axis model is that the early afterglow evolution 
should be charac terized by a rising phase as the jet decelerates and spreads sideways into 



our line-of-sight (IGranot et al. 1120021 ; IWaxmanll2004l ; ISoderberg et al. Il2006el ) 



Another theory suggests that XRFs are intrinsically different from GRBs in their ability 
to couple energy to highly-relativistic material. In this scenario, XRFs are produced in 
explosions characterized b y lower bulk Lorentz factors, 1 < V < 100, than those inferred 
for typical GRBs, Y > 100 ( Zhang. Woosley fc Hegerll2004l ). Thi s may be the result of baryon 
loadi ng of the high-velocity ejecta, a so-called "dirty fireball" (IDermer. Chiang fc Bottcher 
19991 ). Generally speaking, low Lorentz factor explosions may be identified through an 
analysis of their prompt emission sin ce an optically thin sp ectrum at high energies implies 
a lower limit on the Lorentz factor (ILithwick fc Saril 1200 ll ). In the case of X-ray flashes, 
however, there are generally insufficient high energy photons for this type of analysis. For 
these events, detailed modeling of the broadband afterglow may be used to place a lower 
limit on the Lorentz factor. 

Here we present an extensive, multi-frequency data set for XRF 050416a at z = 0.6528 
which extends to t ~ 220 days after the burst. By combining near- infrared, optical, ultra- 
violet, radio and X-ray data we present an in-depth analysis of the afterglow, energetics, 
supernova and host galaxy of XRF 050416a, placing it among the best-studied X-ray flashes 
to date. Moreover, thanks to our dedicated late-time monitoring campaign, we show that 
XRF 050416a produced a strong radio flare at t ~ 40 days accompanied by a brief plateau 
phase in the X-ray band. Finally, we discuss the nature of XRF 050416a in the context of 
off-axis and dirty fireball models for X-ray flashes. 



Observations 



XRF 050416a was discovered by the Swift Burst Alert Telescope (BAT) on 2005 April 
16.4616 UT. The ratio of 15—25 keV and 25—50 keV channel flue nces, j 15-25 kev/Z^s-so kev ~ 
1.1, classifies the event as an X-ray flash ( Sakamoto et al. 2006h. Th is is consistent with the 
low peak photon energy, E p = 15.0 ^2!? keV (jSakamoto et al. 1120061); a 



'actor of ~ 10 lower 



than the typical values observed for long-duration GRBs (jBand et al. II1993T) and a factor 
of ~ 3 larger than the values inferred for XRFs 020903 and 060218 JSakamoto et al. \\200i 
Campana et al. 1120061 ) 



As discussed by ISakamoto et al. I ( 120061 ) , the prompt emission light-curve is character- 
ized by a relatively smooth, triangular peak which is only detected at energies below 50 
keV. The burst duration is T 90 ~ 2.4 sec (15-150 keV), placing it between the classes of 
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short - and long-durat i on G RBs while the hardness ratio shows a clear softening. In addi- 



tion, iSakamoto et al. I (120061 ) note two intriguing features of the data: (1) the rise time of the 
pulse is longer than the decay time, and (2) the cross-correlation lag function (an indication 
of the spectral softening) is —0.0661°;^ sec, apparently inconsistent with the overall soften- 
ing trend observed for the light-curves. This lag function estimate is significantly different 
than the typical values inferred for long-dur ation GRBs and even m ore extreme than the zero 
spectral lags inferred for short-hard bursts (INorris fc Bonnelll 120061 ) . We note, however, that 
the spectral lag for XRF 050416a was derived through a comparison of the two softest BAT 
bands (15-25 and 25-50 keV) and therefore prevents a clear comparison with the spectral lag 
estimates for other Swift GRBs for which 15-25 and 50-100 keV bands are typically used. 
Given that the temporal evolution is not strongly variable, this may indicate that the prompt 
emiss ion was produced by another process (e.g. external shocks, iDermer. Chiang fc Bottcher 



19991 ). Finally, we note that this burst is inconsistent with the lag-luminosity correlation for 



long-duration GRBs which posits that low lum inosity bursts such as XRF 050416a have long 
spectral lags (INorris. Marani fc Bonnelll 120001 ). 



2.1. Early Optical Observations 



Using the roboticized Palomar 60-inch telescope (P60) we initiated observations of the 
field of XRF 050416a at 2005 April 16.4634 UT (2.5 minutes after the burst). In our first 
120 second image we discovered a new source within the Swift/BAT error circle at a = 
12 h 33 m 54 s .58, 5 = +2r03'26".7 (J2000) with an uncertainty of 0.5 arcsec in each coordinate 
based on an astrometric tie to the USNO-B catalog (Figured]). We subsequently monitored 
the afterglow evolution with the Palomar 60-inch, 200-inch, and Siding Springs 2.3-meter 



telescopes in the R- 



and K s — band through t ~ 7 days. 



Aperture photometry was performed on each of the images in the standard method using 
the apphot package within IRAF. Absolu te calibration o f R-, I-, and z' — band light-curves 
was derived using the field calib ration of iHendenl (120051 ) and utilizing the transformation 
equations of ISmith et al. I (120021 ). The -fG-band light-curve was calibrated against 2MASS 
using 15 unconfused sources. The errors resulting from calibration uncertainty (< 10%) 
were added in quadrature to the measurement errors. As shown in Figure [2] and Table [TJ 
the afterglow was I = 18.82 ± 0.11 mag at t ~ 1.6 min (mid-exposure). 

W e supplement these NIR/optical afterglow data with additional measurements from the 
GCNs Li et al. II2OO5I : IPrice et al. II2OO5I: loin et al. ll20od : lToriill2005l : lYanagisawa. Toda fc Kawai 
20051 ) and those reported by lHolland et al. I (120061 ). obtained with the Swift/UVOT and 1.54- 
m Danish Telescope. The resulting dataset spans 1930A (Swift/XJVOT UVW2) to 22000A 
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(K s ); however, we note that the majority of the UV observations are upper limits. Using 
this extended dataset we measure the temporal and spectral properties of the NIR/optical 
afterglow emission. We find the following NIR/optical power-law decay indices (a, where 
F v oc t a ) between ~ 0.01 and 1 day: a Ks = -0.7 ± 0.3, a 7 = -0.7 ± 0.3, a R = -0.5 ± 0.3, 



ay = —1.0 ± 0.3, as = —0.5 ± 0.3, consistent with the values reported by iHolland et al. 
(120061 ). These values imply a mean temporal index of c%iR/opt — —0.7 ± 0.2. 

Finally, we analyze the spectral index (f3 with F v oc v^) within the NIR/optical bands. 
As shown in Figure [3] there are two epochs at which the photometric spectrum is well- 
sampled: t « 0.014 and 0.8 days. We fit each of the observed spectra with a simple power- 
law and find /^NiR/opt ~ —1.3 (xf ~ 0.5) and —1.5 (xf ~ 0.6) for the first and second 
epochs, respectively. As will be discussed in §3.1} the observed steep spectrum is indicative 
of extinction within the host galaxy. 



2.2. Late-time Observations with HST 



Using the Wide-Field Camera (WFC) of the Advanced Camera for Surveys (ACS) on- 
board the Hubble Space Telescope (HST), we imaged the field of XRF 050416a four times, 
spanning 37 to 219 days after the burst. Each epoch consisted of two or four orbits during 
which we imaged the field in filters F775W and/or F850LP, corresponding to SDSS i'- and 
z'- bands, respectively. 



The iKTdata were processed using the multidrizzle routine (IFruchter fc Hookll2002l ) 
within the stsdas package of IRAF. Images were drizzled using pixf rac=0.8 and pixscale=1.0 
resulting in a final pixel scale of 0.05 arcsec/pixel. Drizzled images were then registered to 
the final epoch using the xregister package within IRAF. We astrometrically tied the HST 
and P60 images using 12 unconfused sources in common resulting in a final systematic un- 
certainty of 0.70 arcsec (2a). 

To search for source va riability and remove host galaxy contamination, we used the 
ISIS subtraction routine by lAlardl (120001 ) which accounts for temporal variations in the 
stellar PSF. Adopting the final epoch observations as templates we produced residual images. 
These residual images were examined for positive sources positionally coincident with the 
P60 afterglow. 

Photometry was performed on the residual sources within a 0.5 arcsec aperture. We con- 
verted the photometric measurements to infinite aperture and calculated the corresponding 
AB magnitu des within the native HST filters using the aperture corrections and zero-points 
provided by ISirianni et al. I (120051 ) . Here we made the reasonable assumption that the tran- 
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sient flux is negligible in the template images. For comparison with ground-based data, 
we also converted the F775W measurem ents to Johnson /-ba nd (Vega) magnitudes using 
the transformation coefficients derived by ISirianni et al. I (120051 ) and assuming the F u source 
spectrum implied by the first epoch HST data. 

As shown in Table [2] and Figure HJ the transient is clearly detected in the first epoch 
HST residual images. An astrometric tie between the first and final epochs shows that the 
residual is offset 0.02 ± 0.02 arcsec from the center of the host galaxy. Our residual images 
show the source to be F775W= 24.35±0.02 and F850LP= 23.83±0.03 mag in the AB system 
(/ = 23.82 ± 0.02 mag in the Vega system) at t w 37 days. As shown in Figure [2], these 
values are a factor of ~ 5 above an extrapolation of the early afterglow decay. The observed 
spectral index between the F775W and F850LP filters is Phst ~ —2.8 ± 0.3, significantly 
steeper than the afterglow spectrum observed at early time ( §2.1l and F igure ED as well as the 
typical values measure d for NIR/optical afterglows (PmR/ op t ~ —0.6; iPanaitescu fc Kumar 
20021 ; lYost et al. 1120031 ) . As will be discussed in §4.4[ the timescale and spectral signature of 
the observed flux excess are suggestive of a thermal supernova component. 



2.3. Spectroscopic Observations 



We observed the host galaxy of XRF 050416a with the Low Resolution Imaging Spectro- 
graph on Keck I on 2005 June 6.3 UT (t ~ 50 days). We placed a 1.0 arcsec longslit across 
the host galaxy at a position angle of PA=87°. Data were reduced in standard manner using 
the onedspec and twodspec packages within IRAF. Flux calibration was performed using 
the spectrophotometric standard star BD+284211. 

As shown in Figure [5] and Tabled we detect several strong emission lines in the spectrum 
including H/5, H7, [O II]A3727, and [O III]AA4959, 5006 at a redshift of z = 0.6528 ± 0.0002. 
Adopting the standard cosmological parameters (H = 71 km s _1 Mpc" 1 , Q M = 0.27, Qa = 
0.73), the isotropic gam ma-ray energy release is .E^iso ~ 1.2 ±0.2 x 10 51 erg (1 keV - 10 MeV; 



Sakamoto et al. II2006I ). Compared with typical long-duration bursts, the prompt energy 
release of XRF 050416a is a factor of ~ 100 lower ( Frail et al. 1 2001: B loom. Frail fc Kulkarni 



20031 ; lAmatill2006l and references therein). 
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2.4. Radio Observations 



We began observing XRF 050416a with the Very Large Arraj£] (VLA) on 2005 April 
16.49 UT (t ~ 37 min). No radio source was detected coincident with the optical position 
to a limit of F u < 122 /iJy at v — 8.46 GHz. However, further observations at t w 5.6 days 
revealed a new radio source with F u w 101 ± 34 /iJy coincident with the optical and X-ray 
afterglow positions at a = 12 h 33 m 54.594 s ± 0.002, 5 = 21°03'26.27" ± 0.04 (J2000) which we 
identify as the radio afterglow. 

We continued to monitor the radio afterglow at 1.43, 4.86 and 8.46 GHz through t fu 140 
days (Table HJ). All observations were taken in standard continuum observing mode with a 
bandwidth of 2 x 50 MHz. We used 3C286 (J1331+305) for flux calibration, while phase 
referencing was performed against calibrators J1221+282 and J1224+213. The data were 
reduced using standard packages within the Astronomical Image Processing System (AIPS). 

As shown in Figure the evolution of the radio afterglow is dissimilar from those of 
typical GRBs and inconsistent with a standard blastwave model. Between t ~ 20 and 40 
days an abrupt rebrightening (factor of ~ 3) is observed at all radio bands with a temporal 
index steeper than a ra d ~ 0.9. Following this radio flare the emission decays rapidly with an 
index of a ra d < —1.5. The peak radio l uminosity at 8.46 GHz, L rad 1.3 x 10 31 erg s" 1 Hz" 1 , 
is typical for GRBs (IFrail et al. 1120031). a factor of ~ 10 and 10 3 larger than those of XRFs 



020903 and 060218, respectively flSoderberg et al. ll2004aL l2006d ). and between 10 2 and 10 6 
times hig her than the pea k radio luminosities observed for optically-selected Type Ibc su- 
pernovae (ISoderbergi 120061 ) . By t ~ 105 days the radio afterglow is no longer detected at 
any frequency. 



2.5. X-ray Observations 



The afterglow of XRF 050416a was observed with the Swift X-ray Telescope (XRT) be- 
ginning 1.2 minutes after the burst an d continuing through t ~ 69 days, placing it among the 



best studied X-ray aftergl ows to date (INousek et al. M2006I : lO'Brien et al. II2006I ). From their 



analysis of the XRT data, iMangano et al. I ( 2006 ) report a spectral index of f3 



6.8 



+1.0 
--L2- 



X 10 



cm 



.1 04+ - 11 
corre- 



and evidence for significant absorption in the host galaxy, iV# 
sponding to E(B — V) TCSt = 1.2 ± 0.2 assuming the conversion of iPredehl fc Schmittl (jl995l ) 



1 The Very Large Array and Very Long Baseline Array are operated by the National Radio Astronomy Ob- 
servatory, a facility of the National Science Foundation operated under cooperative agreement by Associated 
Universities, Inc. 
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and a Milky Way extinction curve (|Peilll992l ). As shown in Figure |7J at t — 10 hours the 
X-ray luminosity was Lx,i SO (t = 10 hr s) ^2.3 x 10 45 erg s" 1 , placing it at the lower edge of 
the o bserved distribution for GRBs (IBerger. Kulkarni &: Frail! 120031 ; iFreedman fc Waxman 
200lh . 



As discussed by iMangano et al. I (120061 ) and shown in Figure [7J the early evolution of 
the X-ray afterglow can be characterized by three phases: (1) an initial steep decay (2) a 
flattening between t w 7 and 20 minutes during which the flux is roughly constant, and (3) 
a resumed decay through t ~ 20 days. These three p hases have been shown to be ubiquitous 



among Swift X-ray afterglows (INousek et al. II2006I ). However, at t ~ 20 days there is a 
second flattening which continues through t ~ 40 days. By 70 days the X-ray afterglow 
is no longer detected, implying a significant steepening to ax ^ — 1.8 between the last two 
observations. We note that the timescale for the X-ray flattening and subsequent steep decay 
is coincident with the observed radio flare. 



3. Properties of the Early Afterglow 



Using the detailed multi-frequency observations of the XRF 050416a afterglow we can 
constrain the physical properties of the ejecta and the circumburst density. We adopt a 
standard relativistic blastwave model in which the afterglow emission is produced through 
the dynamical interaction of the ejecta with the surrounding medium (the forward shock, 
FS) with an additional component from shock heating of the ejecta (the reverse shock, RS). 
In this scenario, the total energy density is partitioned between relativistic electrons, e e , 
and magnetic fields, e r , while the thermal energy of the shocked protons accounts for the 
fraction remaining (see |Piranlll999l for a review). The shocked electrons are accelerated into 
a power-law distribution, N(j) oc 7 _p above a minimum Lorentz factor, 7 m . The emission 
resulting from the forward and reverse shock components is described by a synchrotron 
spectrum characterized by three break frequencies — the self-absorption frequency, u a , the 
characteristic frequency , u m , and the cooling frequency, u c — and a flux normalization, F„ m 
( iSari. Piran fc Narayanl I1998T). In modeling the afterglow spectral and temporal evolution, 
we adopt the formalism of iGranot fc Saril (120021 ) for a relativistic forward shock expanding 
into a constant density circumburst medium. 
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3.1. Preliminary Constraints 

In fitting the forward shock model to the afterglow data of XRF 050416a we use only 
observations between 0.014 and 20 days when the afterglow follows a simple power-law 
evolution. To constrain the spectrum of the forward shock, we first investigate the afterglow 
evolution in the optical and X-ray bands. As shown in Figure t he X-ray data b e tween 



0.014 and 20 days are reasonably fit with ax ~ —1.1 (xf ~ 0.70). iMangano et al. I (120061 ) 
report that the X-ray spectral index on this same timescale is (3x ~ —1.04, leading to 
a — 3/5/2 « 0.5. A comparison to the standard closure relations, a — 3/3/2 = {y m < v < u c ) 
and a — 3(3/2 = 1/2 [y > u c ) indicates that v x > v c . This conclusion is supported by 
the near-IR to X-ray spectral slope, (3k,x — —0.47 ± 0.06 at t ~ 0.6 days, which is flatter 
than (3x as expected if ^NiR/opt < v c < Vx- Therefore, the X-ray observations suggest that 
p = -2/5 » 2.1. 

Next we consider the spectral index within the NIR/optical bands. As discussed in 
§2.11 and shown in Figure [31 the observed NIR/optical spectral index on this timescale is 
AvriR/opt ~ —1.3 to —1.5. These values are significantly steeper than (3k,x and imply that 
the optical flux is suppressed by host galaxy extinction. Making the reasonable assumption 
that v m < ^NiR/opt on the timescale of our afterglow observations, and adopting p ~ 2.1 
as indicated by the X-ray observations, we estimate /?NiR/ op t = —(p — l)/2 ~ —0.55 for the 
intrinsic spectral index of the NIR/optical afterglow. Adopting this value for /?NiR/opt, we find 
that 



joth NIR/optical spectra are reaso nably fit with a Galactic extinction of E(B — V) 



0.03 flSchlegel. Finkbeiner fc Davislll998l ) and a host galaxy component of E(B—V) 



est 



0.28 



(Figure [3]). Here we have assumed a Milky Way extinction model for the host (|Peilll992l ) 
but note that a comparable fit may be obtained for an SMC extinction curve. We further 
note that this optically-derived extin ction estimate is lower than that inferred from the X- 



GRBs iGalama & Wiierd (J200l|). 



ray sp ectrum (IMangano et a/. 1 120061 ) . consistent with the trend observed for long-duration 



With this extinction correction, the near-IR to X-ray spectral index becomes (3k,x ~ 
— 0.5±0.1, consistent with our estimate for the intrinsic spectral index within the NIR/optical 
band. Moreover, the extinction-corrected NIR/optical spectral index and observed average 
temporal index of aNiR/ op t = —0.7 ± 0.2 are consistent with the standard closure relation: 
a — 3/3/2 = ~ 0.1. We also note that this supports our assumption of a constant density 
medium since in a wind environment the expected temporal index is steeper than a = —1.25, 
and thus inconsistent with the observed values. Using all the available optical and X-ray 
observations we estimate that v c » 1 x 10 17 Hz at t = 1 day. 

Next we compare the near-IR and radio afterglow data to constrain v m and the peak 
spectral flux, F Um . Assuming that v m passed through the NIR/optical bands near the time of 
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our first i^ s -band observations implies that the peak spectral flux is roughly comparable to 
the extinction-corrected K s flux: Fk s ~ 230 /xJy at t ~ 11 minutes. Here we focus on the K s - 
band data since they are the least sensitive to host galaxy extinction, which we estimate to 
be A K « 0.24 mag (a 20% increase in flux) for E(B — V) Icst ~ 0.28. Scaling these constraints 
to t = 1 day [y m oc t" 1 ' 5 and F Vm oc t°) and accounting for the smooth shape of the spectral 
peak, we find u m 4.0 x 10 11 Hz and F Vm ~ 350 [iJy. Here and throughout, F Vm is the 
asymptotic extrapolation of the smooth spectrum peak and is therefore slightly higher than 
the intrinsic peak flux. We note that since the NIR/optical data require F Vm oc v~(p~ 1 )/ 2 , 
lower values of v m imply increasingly higher values of F Vm at the time of the first f^-band 
observations. 

Finally, we test that these constraints are consistent with the radio observations. Given 
the evolution of u m , these constraints predict that the spectral peak should pass through the 
radio band at t ~ 13 days with an extrapolated peak flux density of F Vm ~ 350 /iJy, roughly 
consistent with the 4.86 GHz observations on this timescale. We emphasize that the early 
steady decay of the NIR/optical data require that v m passes through the radio no later than 
13 days. Finally, we note that the radio spectrum is optically thin throughout the timescale 
of VLA monitoring (see Figure ED and thus we observationally constrain v a to be below 1.43 
GHz. 



3.2. Forward Shock Broadband Model 



Adopting these constraints we apply a broadband afterglow model fit to the multi- 
frequency data in order to determine the physical parameters of the burst. The four spectral 
parameters (F Vm , u a , v m and v c ) are fully determined by four physical parameters: the 
isotropic ejecta energy, Ek,sso, the energy density partition fractions, e e and €b, and the 
circumburst density, n. Therefore by constraining the four spectral parameters through 
broad-band obs ervations, we are able to de ter mine a uniqu e solution for the four physical 
parameters (see lSari. Piran &: Narayanlll998l and lPiranlll999l for reviews). Although the radio 
observations provide only an upper limit on u a , we are able to define a range of reasonable 
values by requiring that e e ,es < 1/3 which accounts for an equal or greater contribution from 
shocked protons. This requirement excludes unphysical solutions in which the sum of the 
contributions from shocked electrons, protons and magnetic fields exceed the total energy 
density. Combined with the observed constraints for F Vm , v m and v c we find the following 
ranges for the physical parameters: 



E KM « (8.2 - 14) x 10 51 erg 
n w (0.33 - 4.2) x 10" 3 cm" 3 



(la) 
(lb) 



e e » (0.20 - 1/3) 
e B « (0.072 - 1/3). 



(lc) 
(Id) 



As shown in Figures [21 [BJ and [7J this model provides an adequate fit to the broadband data 
between t ~ 0.01 and 20 days. 



3.3. Collimation of the Ejecta and Viewing Angle 



The lack of an observed jet break in the X-rays to t ~ 20 da ys, together with the inferred 



to j « 3.1tf 8 ^ is0j52 n_ 



Sari. Piran fc Halpernlll999l ) 



physical parameters constrain the opening angle of the jet (e.g 

1 N ""'Yg (1 + z)~ 3 ^ 8 > 6.9 degrees. Here, tj is the jet break time in days 



of the jet opening angles inferred for lone-duration GRBs, Q~ ~ 5° ( 


31oom. Frail & Kulkarni 


2003; 


Ghirlanda, Ghisellini & Lazzati 


2004; 


Soderbere et al. 


2006a 


and references therein). 



This indicates that the beaming-corrected ejecta energy release is E K = E Kiso (l — cosO) ~ 
(9.8 x 10 49 — 1.4 x 10 52 ) erg where the range includes the uncertainty in Ek^ and the lower 
limit on 



j. Moreover, we expect the blastwave to become n on-relativistic on a timescale 

V3. 



2.0 E 5i -n_ 3 



1/3 



(0.6 — 6.8) yrs flLivio fc Waxmanll2000l ). On a similar timescale, the 



ejecta are predicted to approach spherical symmetry after which the blastwave evolution 
is we l l described by the Sedov-yon Ne umann- Taylor (SNT) solution (jZel'dovich fc Raizer 
20021 ; iFrail. Waxman fc Kulkarnilboooh ; in this regime the afterglow emission decays with 
-9/10 ( — 1) for frequencies below (above) the cooling frequency. 



a 



Next, the early steady decay of the X-ray and NIR/optical afterglow indicates that 
the jet collimation axis is directed roughly along our line-of-sight. In comparison, GRBs 
viewed significantly off- axis (6^ > 20 j) are predicted to show a rising afterglow l i ght-curve 



as the jet spreads sideways an d intersects our viewing angle ( IGranot et al. 



2002: Waxman 



20041 ; ISoderberg et al. Il2006el ). Here, # b s is the angle between our line-of-sight and the jet 
collimation axis. We conclude that the ejecta are viewed roughly on-axis and therefore the 
inferred beaming-corrected energies are not affected significantly (if at all) by viewing angle 
effects. 



4. Properties of the Late-time Afterglow 

Next we address the nature of the late-time broadband afterglow evolution with spe- 
cial attention to the strong radio flare observed at ~ 40 days. Radio flares have been 
noted for several other GRBs, though only at early times (e.g. GRB 990123 at t < 1 day; 
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Kulkarni et al. Ill999l ). Based on their observed timesc ale and evolution, r adio flares are typ- 



ically attributed to emission from the reverse shock (ISari &: Piranlll999l ). Here we present 
detailed radio observations for XRF 050416a which show for the first time a strong radio 
flare at late time. Possible causes for a late-time radio rebrightening in clude the emis- 



sion from a dec elerating jet initially directed aw ay from our l ine-of-sight rtWaxman 



Li fc Song! 12004 ) , circumburst density variations (IWijersI l200ll ; iRamirez-Ruiz et al. 



2004 



200lh . 



and energy injection fr om a slow shell catching up to and refreshing the afterglow shock 
( iRees fc Meszaroslll998l ). We discuss each of these possibilities below. 



4.1. Off-axis Jet Emission 



It has been shown that the observational signature of a relativistic jet viewed signifi- 
cantly away from our line-of-sight is a rapid achromat ic rise in the early afterglow light-curves 
( ]Paczynskill200ll ; iGranot et al. 1120021 ; IWaxmanl I2004T ) . In this scenario, the observed peak of 
the afterglow emission occurs as the spreading jet crosses our viewing angle. The timescale for 
this peak is ~ 100 days for a GRB jet with typical paramete rs (Ek = 10 51 erg , n = 1 cm -3 , 
e e = €b = 0.1, 9j = 5°) viewed from an angle ^obs = 30° ( ISoderberg et al. Il2006el ). The 



subsequent afterglow evolution is the same as that seen by an on- axis observer, decayin g 
steeply with a = —p ~ [—2 to —3] for frequencies above v m ( ISari. Piran fc Halpernlll999l ). 
The observed timescale and evolution of the XRF 050416a radio flare are therefore roughly 
consistent with the predictions of an off-axis relativistic jet. 



However, as noted by lMangano et al. I (120061 ) and discussed in §3.31 the early and steady 
decay of the XRF 050416a broadband afterglow implies the presence of relativistic ejecta 
directed along our line-of-sight. These ejecta are also responsible for the production of the 
observed prompt emission. Attributing the strong late-time radio flare to an energetic off-axis 
relativistic jet therefore implies that multiple relativistic ejecta components were produced 
in the explosion. Moreover, the steep rise and peak flux of the radio flare imply sharp edges 
for the off-axis jet and a kinetic energy 2 to 3 times larger than that of the on-axis ejecta. 
While this scenario cannot be ruled out, we consider it unlikely given the complicated ejecta 
geometry required. 



4- 1.1. A Receding Jet 



Li &: Songi (120041 ) describe a related scenario in which a strong late-time radio flare 
is observed from a receding jet initially directed anti-parallel to our line-of-sight. In this 
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case, early afterglow emission is expected from the approaching jet while emission from the 
receding jet is expected on a timescale 5 x t NR due to the light travel time delay. This 
scenario is appealing in that it may explain both the early and late-time afterglow emission 
observed for XR F 050416a w ithin the standard framework of engine-driven (accretion-fed 
compact source; |Piranlll999l ) double-sided jets. However, as discussed in §3.31 broadband 
modeling of the afterglow emission predicts the non-relativistic transition to occur no earlier 
than 0.6 yrs - too late to explain the radio flare at 40 days. Moreover, the peak emission 
from a receding jet should be comparable to the r adio flux at £nr arid decay with a temporal 
index given by the Sedov solution, a ra d = —9/10 (ILi fc Songj|2004l ); both of these predictions 
are inconsistent with the observations. We conclude that afterglow emission from a receding 
jet is unlikely to produce the observed radio flare. 



4.2. Circumburst Density Variations 



It has been argued that abrupt variations in the circumburst medium can produce a 
strong rebrightening in the radio afterglows of GRBs. Specifically, the dynamical inter- 
action of the forward shock with a wind-termination shock at r ~ 1 pc is predicted to 



cause a rebrightening of t he ra dio afterglow on a timescale of a few years (jWijerd 12001 



Ramirez-Ruiz et al. I l200ll . 120051 ) . Here it is assumed that the blastwave is expanding non- 

relativistically when it encounters the d ensity jump, overall consistent with the observed 

timescales for non-relat ivistic transitions ( Frail. Waxman fc KulkarnlboOO ; Berger. Kulkarni fc Frail 
20041 ; iFrail et al. II2005I ). For comparison, the interaction of a relativistic blastwave with an 
abrupt (step function) density jump is not expected to cause strong variations in the af- 
terglow ligh t-curves, where here it has been assumed that the post-jump expansion is also 
relativistic (INakar fc Granotll2006l ). 



In the case of XRF 050416a, the strong radio flare and X-ray plateau phase occur on a 
timescale when the blastwave is still relativistic, and therefore an abrupt circumburst den- 
sity jump appears an unlikely explanation. This is supported by the fact that during the 
relativistic regime, the flux at frequen cies above v c should be unaffected by circumburst 
density variations (IGranot fc Sar We speculate, however, that a very large density 

jump may be able to decelerate the blastwave to non-relativistic speeds on a very short 
timescale and may therefore be able to explain the unusual late-time afterglow evolution. 
The increase in density would cause a shift in v a which may explain the peculiar evolution 
of the spectral index in the radio band (Figure [6]). Based on our afterglow modeling ( §3. 2[) . 
we estimate that the circumburst radius of the forward shock at the onset of the radio 
flare was roughly ~ 1 pc, roughly consistent with the radius of a wind-termination shock 
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( IGarcia-Segura. Langer fc Mac Lowlll996l ; IChevalier. Li fc Franssonll2004l ). In this scenario, 
we expect a post-jump self-similar evolution consistent with the Sedov- Taylor solution. Since 
the late-time radio and X-ray data are not sufficiently sensitive to trace the post-flare evo- 
lution, we cannot rule out this possibility. 

Finally, we investigated a scenario where the radio/X-ray flare is produced by the dy- 
namical interaction of the quasi-spherical, non-relativistic SN ejecta with a CSM density en- 
hancement. In fact, density jumps have been invoked to explain radio modula tions observed 
for local (non-relativistic) SNe Ibc (e.g. SN2003bg; ISoderberg et al. Il2006bl ). Adopting a 
simple minimum energy calculation ( Kulkarni et al. 19981 ) and requiring that the shock en- 
ergy is equally partitioned between magnetic fields and relativistic electrons, we find that 
attributing the strong radio emission to the quasi-spherical SN component requires that the 
SN ejecta is relativistic (T ~ 10) at the time of the radio/X-ray flare. However, as discussed 
above, strong flux variations are not expected while the blastwave is relativistic. Combined 
with the fact that the flare is at least a factor of 10 2 more radio lumin ous than any other 
SN Ibc ever observed (including GRB-SN 1998bw. Kulkarni et al. II1998I ). we conclude that 
the radio/X-ray flare can not be attributed to the quasi-spherical SN ejecta. 



4.3. Energy Injection 



An episode of energy injection may also cause a rebrightening of the afterglow flux. 
Energy injection may arise from long-lived central engine activity or under the framework 
of a "refreshed" shock where a slow moving shell ejected du r ing the initial burst eventuall y 
catches up with the afterglow shock (IRees fc Meszaroslll998l ; iGranot. Nakar fc Piranll2003l ). 
The observed ratio of the settling time to the epoch of injection, 5t/t, can distinguish between 
these two scenarios. For 5t/t < 1 the injection is produced by engine activity, while 5t/t > 1 
indicates a refreshed shock. In the case of XRF 050416a, the X-ray and radio data after 20 
days imply 5t/t > 1, suggesting that the ejecta were refreshed on this timescale. 

As shown in Figure El an extrapolation of the early radio evolution lies a factor of ~ 3 
below the observed flux at the onset of the radio flare. An energy increase affects the spectral 
parameters according to the following scalings: v a oc E 



and v c oc E K . so 



1/5 
Kasoi P m 



r,l/2 
v if,iso> 



oc E 



1/2 

(IGranot fc Saril 120021 ) . Given that v m is within the radio band on this 
timescale we have F^ ra d oc Ek^o and thus the radio flare corresponds to an energy injection 
of a comparable factor, ~ 3. For the X-ray band, F Vy x oc E^^J^, thus for p ps 2.1 an 
energy injection of a factor of ~ 3 corresponds to a comparable increase in the X-ray flux. 
The Swift/XKT observations suggest a flattening on this timescale. Here we adopt the 
conservative assumption that the shock microphysics (e e , e# and p) do not evolve during the 
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energy injection. 



One important implication of the energy injection model is that the post-injection 
asymptotic temporal decay should be the same as that before the injection. This predic- 
tion is consistent with the energy injection episodes invoked for GRBs 02100 4 and 030329 
( iBjornsson. Gudmundsson fc Johannessonl 12004 ; iGranot. Nakar fc Piranll2003l ). However, in 
the case of XRF 050416a, the late-time radio and X-ray data suggest a steep post-injection 
decay, a > —2, significantly steeper than that observed pre-injection. Moreover, the af- 
terglow should asymptotically approach a flux normalization larger by a factor of ~ 3 in 
both the radio and X-ray bands. As shown in Figures [6] and [7J this appears inconsistent 
with the observations which suggest that the late-time afterglow is comparable (or fainter) 
than an extrapolation of the early afterglow model. We note, however, that the faintness of 
the late-time afterglow and the steep post-injection decay may be explained if a jet break 
occurred on roughly the same timescale as the energy injection. 



In comparison with other late-time afterglow studies (IFrail. Waxman &: K ulkarni 



2000 



Tiengo et auhooA berger. Kulkarni fc Frai3l2004J : iKouveliotou et al. 1120041 : IFrail et al. Il2005h , 
the radio flare and X-ray flattening observed for XRF 050416a are clearly atypical for long- 
duration GRBs. We therefore attribute the observed evolution to an unusual scenario in- 
volving either a large circumburst density jump or a late-time injection of energy (from a 
slow shell or off-axis ejecta). 



4.4. An Associated Supernova 

The HST measurements at late-time can be used to constrain any possible contribution 
from an associated supernova. Based on previous studies of GRB-SNe, the thermal emission 
from an associated supernova is predicted to reach maximum li ght on a timescale of 20x (l+z) 



days with a peak magni tude of My rest fa —20 mag or fainter (IZeh. Klose fc Hartmannl 12004 



Soderberg et al. I l2006dl and references therein). Observationally, the emergence of a SN 
component produces a steepening of the optical spectrum as the supernova nears maximum 
light and dominates the afterglow emission. 

An extrapolation of the broadband afterglow model to 40 days shows that the HST 
data are brighter by a factor of ~ 5. For comparison, the radio flare and X-ray plateau on 
this same timescale represent flux density enhancements by factors of roughly ~ 3 and 2, 
respectively. 

As discussed in §2.21 we measure the spectrum of the optical transient within the ACS 
bands and find (3hst ~ —2.8 ± 0.3 at 37 days. After correction for extinction (see §2.11) . the 
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implied spectral index becomes (3hst ~ —1-9. For comparison, the spectral index between 
the radio and X-ray bands on this timescale (coincident with the radio flare and X-ray 
flattening) is /3 ra ,dx ~ —0.56 ± 0.04 and the indices w ithin the bands are /3 ra d ~ —0.6 ± 
0.2 (Figure EJ) and (3x ~ —1-0 (IMangano et al. II2006I ). respectively. Given that the HST 
derived spectral index, even after correction for extinction, is (1) significantly steeper than 
the NIR/optical index observed at early time, /^NiR/opt ~ — 1-4, (2) significantly steeper than 
A-adx measured at a comparable epoch, and (3) inconsistent with t he range of synchrotron 
spect ral indices predicted for a relativistic blastwave (J3 = [—1.5, 2.5]. ISari. Piran fc Narayan 
19981 ). we conclude that the optical flux at t ~ 40 days is dominated by another emission 
process, likely an associated SN. 

Next, to determine if the late-time optical data are consistent with the temporal evolu- 
tion of a typical GRB/XRF-associated supernova, we compare the HST flux values with syn- 
thesiz ed SN light-curves. We compiled UBVRIJHK observations of SN 199 8bw from the liter- 
ature (IGalama et al. Hl998l ; lMcKenzie fc Schaeferlll999uPatat et al. 1120011) and smoothed the 
extin ction-corrected (Galactic component of E(B-V)=0.059, ISchlegel. Finkbeiner fc Davis 
19981 ) light-curves. Here we ass ume that the SN 19 98bw host-galaxy extinction is negligible, 
consistent with the findings of iPatat et al. I ( 120011 ) based on a spectroscopic analysis of the 
host galaxy. We then produced fc-corrected NIR/optical light-curves of SN1998bw at the 
redshift of XRF 050416a by interpolating over the photometric spectrum and stretching the 
arrival time of the photons by a factor of (1 + z). 

Shown in Figure[2]are the synthesized light-curves for SN 1998bw at z — 0.6528, summed 
together with the afterglow model. The HST data are roughly comparable with the flux 
normalization and evolution of the summed model. Therefore, the temporal and spectral 
properties of the HST data suggest that XRF 050416a was associated with a supernova 
similar to SN 1998bw. However, we caution that the temporal coincidence of the SN peak 
with the radio flare makes it difficult to estimate the relative contributions of the SN and 
afterglow. 



5. Host Galaxy Properties 



We now turn to the properties of the GRB host galaxy. We measure the brightness of 
the host galaxy in the final HST epoch to be F775W= 23.1 ±0.1 mag (I = 22.7 ± 0.1 mag). 
These values are not corrected for extinction. At z — 0.6528 the rest-frame 5-band is traced 
by the observed F775W band, leading to an absolute magnitude, Mb ~ —20.3 ± 0.1 mag, 
or a luminosity Lb » 0.5L*. T his host luminosity is similar to t hat inferred for XRF 03 0528 
flRau. Salvato fc Greinerl l2005h and the hosts of typical GRBs Le Floc'h et al. Il2003h . At 
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z = 0.6528, the measured offset of the optical transient ( §2.ip relative to the center of the 
host galaxy corresponds to 140 ± 140 pc. This offset is a factor of ~ 10 s maller than the 
median value for long-duration GRBs (IBloom. Kulkarni &: Djorgovskill2002l ). 



As shown in Figure [5j the host exhibits several emission lines typical of star-forming 
galaxies. We estimate the star formation rate in the host galaxy from the observed fluxes 
of the various emission lines. Using the flux of the [O II|A3727 l ine, F[on] ~ 9.6 x 10~ 17 
erg cm" 2 s" 1 (Table [3J, and the conversion of IKennicuttl (119981 ). SFR = (1.4 ± 0.4) x 
10 _41 L[oii] M yr _1 , we find a star formation rate of about 2.5 ± 0.7 M Q yr _1 . From the 
H/3 line flux, F^p ~ 3.7 x 10 -17 erg cm -2 s" 1 , and assuming th e Case-B recombination 
ratio of Fn a / F^p = 2.87 and the conversion of IKennicuttl (119981 ). we infer a star forma- 
tion rate of SFR = 7.9 x 10~ 42 L Ha « 1.5 ± 0.2 M yr" 1 . Thus, we conclude that the 
star formation rate (not corrected for extinction) is roughly 2 M Q yr _1 . We note that 
the observed ratio of H/y/H/S = 0.3 ± 0.1 compared to the theoretical value of about 0.47 
(IOsterbrocklll989l) suggests a significant extinction corr ection (factor of ~ 10) following the 



method of (jCalzetti. Kinney fc Storchi-Bergmannlll994l ). We conclude that the star forma- 
tion rate for the host of XRF 050416a is similar to those in ferred for long-duration GRBs 



(jChristensen. Hjorth fc Gorosabelll2004l ; iBerger et al. 



tude larger than that inferred for XRF 020903 (jSoderberg et al. Il2004al ; iBersier et al. II2006I ). 



2003aT) and at least an order of magn i 



The combination of the inferred star formation rate and host luminosity provides a 
measure of the specific star formation rate. We find a value of 4 M yr _1 L*/L (uncorrected 
for extinction), which is about a factor of two lower than the mean specific star formation 
rate for the hosts of long-duration GRBs (jChristensen. Hjorth fc Gorosabell 12004 ) . 



Next, we use the relative strengths of the oxygen and hydrogen emission lines to infer 
the ionization state and oxygen abundance. The relevant indicators are R23 = l og (-Fjou] + 
F[cnr(\ / Fug) ~ 0-74 and = log ( - ^[ottt] / -^[qtt]) ~ 0.055. Using the calibrations of iMcGaugh 



(Il99ll ) and iKobulnicky fc Kewleyl (120041 ) we find that for the upper branch the metallicity 
is 12 + log(0/H) m 8.6 while for the lower branch it is about 7.9; the two branches are 
due to the double- valued nature of R23 in terms of metallicity. Thus, t he host metallicity of 
XRF 050416a is 0.2 - 0.8 Z Q , larger than that inferred for XRF 060 218 JModiaz et al 



2006; 



Wiersema et al. 1 120071 ) and comparable to that for XRF 030528 (IRau. Salvato &: Greiner 
20051 ) . Moreover, this range is somewhat higher than the typic al metallicities for GR B 
hosts, some of which have metallicities that are ~ 1/10 solar (e.g.. iProchaska et al. II2006I ). 
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Is XRF 050416a a Short Burst? 



With a prompt emission duration of just T 90 ~ 2.4 sec ([Sakamoto et al. 1120061 ). it is 
interesting to consider XRF 050416a as a member of the short-hard class of gamma-ray 
burst s, popularly believed to result from the coalescence of neutron stars or black holes 



|e.g., lEichler et al. Ill989l ). Based on the bi- modal BATSE duration distribution, bursts 
with T90 > 2 seconds are assumed to belong to the long-duration class (jKouveliotou et al. 



19931 ). although a decomposition of the overlapping distributions suggests that a small frac- 
tion of short-hard bursts (SHBs) have durations longer than this cut-off. The distinction 
between short-hard and long-duration bursts is further complicated by the detection of soft 



X-ray tails lasting sever al seconds following SHBs 050709 and 050724 (jVillasenor et al. 112005 



Barthelmy et al. II2005I ). This suggests th at SHBs are not necessarily characterized by a pure 
hard emission spectrum (ISakamotoll2006l ) . Related to this is sue is the use of spectral lags to 
distinguish between long and short bursts. As discussed by iNorris fc Bonnelll (120061 ). long- 
duration bursts typicall y have longer lags that correlat e with isotropic equivalent prompt 
gamma-ray luminosity (INorris. Marani &: Bonnelll |2000| ) . On the other hand, SHBs have 
negligible (or even negative) spectral lags. 

In the case of XRF 050416a, the prompt duration places it between the long- and short- 
duration classes. The low value of E p may suggest that it belongs to the long-duration class, 
however, it is beco ming clear th at hardness cannot be used to reliably distinguish between 



the two classes (see iNakarl 120071 ). Similarly, while the negative lag inferred for XRF 050416a 



may suggest a SHB classification, ex amples do exist of long-duration BATSE bursts with 
negative lags (INorris fc Bonnelll 120061 ). 



Afterglow modeling may provide additional clues. The range of values we infer for the 
XRF 050416a beaming-corrected energies are overall consistent with those of long-duration 
GRBs including the subclass of sub-ener getic bursts. At the same time, they are roughly con- 



sistent with the values inferred for SHBs (IFox et al. ll2005l : lBerger et al. ll2005l : lSoderberg et al. 
2006al ). However, the low circumburst density, n ~ 10~ 3 , is a factor of 10 to 10 4 smaller 



Yost et al. 



than the typical values inferred for long-durati on GRBs (e.g. iPanaitescu fc Kumai 



2005 



2003: 



Cheva 



Berger et al. 



2005 



ier. Li &: Fransson 



Bloom et al. 112006b 



2004) but comparable to those of SHBs (IFox et al 



Burrows et al. 1120061 ; |Panaitescull2006l ; ISoderberg et al. 



2002; 



2006al ). The observed radio flare and X-ray flattening at 20 days are atypical for both long- 
and short- duration GRBs and therefore cannot be used to classify this event. However, 
it is interesting to note that large variations in the circumburst density are more natu- 
rally explained in the context of a massive stellar progenitor with in teracting stellar winds 



( jGarcia-Segura. Langer fc Mac Lowill996l ; iRamirez-Ruiz et al. Il200ll ). 

Next is a discussion of the XRF 050416a host galaxy since long- and short- bursts 
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may also be distinguished by their environments (see iBergerl 120061 for a recent review). 
As discussed in $51 the host is a star-forming galaxy with an inferred SFR and metallic- 
ity comparable to those of long-duration GRBs. Moreover, XRF 050416a is located near 
the center and brightest part of its host galaxy. It is thus consistent with the locations 



of long-duration bur sts with respect to their hosts (IBloom. Kulkarni fc Djorgovskil 12002 



Fruchter et al. 1 120061 ) . For comparison, SHBs are typically localized to low SFR hosts with 
significant old stellar populatio ns at radial offsets up to a factor of 10 larger than those 



of typical long-duration bursts (IBerger et al. 
Soderberg et al. Il2006al ). 



20051 : iFox et al. 1 120051 : IBloom et al. I l2006al : 



Finally, the discovery of an associated Type Ic supernova is perhaps one of the best meth- 
ods to distinguish between long- and short-bursts. The discovery of several long-duration 
bursts and X-ray flashes at z < 0.3 in the last few year s has firmly established that GRBs 



and XRFs are accompanied by supernovae of Type Ic (IGalama et al. 1 1 19981 ; IStanek et al. 



2003: Malesani et al. 



2004 ; ISoderberg et al. Il2005bl ; IPian et al. 1120061 ) . A study of the lumi 



nosity distribution for GRB/XRF-SNe reveals a significant d ispersion, implying a spread of 
(at least) an order of m agnitude in peak optical luminosity ( IZeh. Klose &: Hartmannl I20Q4J : 
Soderberg et al. Il2006dl ). At the same time, deep imaging of SHBs has co nstrained any as 
sociated SN emission to be up to ~ 100 times less luminous as SN 1998bw (IFox et al. 112005 



Soderberg et al. Il2006al ). In the case of XRF 050416a, the temporal and spectral properties 
of the HST data suggest that XRF 050416a was accompanied by a supernova with a peak 
luminosity roughly similar to SN 1998bw. We stress, however, that the temporal coincidence 
of the radio flare with the SN peak complicates any study of the SN properties. 

Based on the associated SN and large-scale environmental properties, we conclude that 
XRF 050416a is a member of the long-duration class of gamma-ray bursts. This event high- 
lights the difficulty in classifying bursts based on their prompt emission properties alone. 



7. Discussion and Conclusions 



We present extensive broad-band NIR/optical and radio data for the afterglow of the 
X-ray flash 050416a and show that it is localized to a star-forming (> 2 M yr _1 ) host galaxy 
at z = 0.6528. Along with XRFs 020903, 030723 and 060218, this burst is one of the best 
studied X-ray flashes to date. Moreover, XRF 050416a is only the third XRF with a spectro- 
scopic redshift for which a broadband afterglow study has been performed and the physical 
parameters have been constrai ned. The isotropic-equi valent prompt and kinetic energy re- 
leases are -E 7i i SO ~ 1.2 x 10 51 ( jSakamoto et al. 1 120061 ) and Ek,\ so ~ 10 52 erg respectively. 



These values are a factor of 10 2 times larger than those of XRF 020903 (jSakamoto et al. 
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2004 ; ISoderberg et al. I l2004al ) and up to 10 4 ti mes larger than those inferred for sub- 



energetic GRBs 980425, 031203 and XRF 060218 ( 


Pian et al. 


2000l: 


Kulkarni et al. 


1998 




Sazonov. Lutovinov & Sunvaev 


2004 


Soderberg et al. 


2004b 


Campana et al. 


2006 


Soderberg et al. 



2006dJ. 



Adopting the results of our standard synchrotron model and tj > 20 days we constrain 
the collimation of the ejecta to 9j > (6.9 — 9.7)°, slightly larger than the median value o f 5° 
for typical GRBs (IBloom. Frail fc Kulkarni! 120031 ; iGhirlanda. Ghisellini fc Lazzatill2004l and 
references therein). This indicates that the beaming-corrected energy release is Ek ~ (9.8 x 
10 49 — 1.4 x 10 52 ) erg and & (8.6 x 10 48 — 1.2 x 10 51 ) erg, implying a total relativistic energy 
yield of E tot (1.1 x 10 50 — 1.5 x 10 52 ) erg, which straddles the median value for cosmological 



GRBs, E tot m 2 x 10 51 erg (IBerger et al. Il2003bl ). However, the efficiency in converting the 
energy in the ejecta into 7-ra ys is 77 7 = E^/(Ek + E 7 ) (0.04 — 0.08), signifi cantly lower than 



the typical values for GRBs (IPanaitescu fc Kumarl 120 02; Yo st et al. 1120031 ) and comparable 



to that inferred for XRF 020903 (r? 7 « 0.03. ISoderberg et al. Il2004af ). This strengthens the 



idea that XRFs and GRBs are distinguished by the i r ability to couple significant e nergy to 
highly-relativistic material ( ISoderberg et al. Il2004al : IZhang. Woosley fc Hegerll2004j ). 



In addition to the burst energetics, several key results emerge from our broadband 
analysis of XRF 050416a. First is the detection of a bright, late-time radio flare accompanied 
by an observed flattening in the X-ray bands which we attribute to a large circumburst 
density enhancement or episode of energy injection (refreshing shell or off-axis ejecta) at t ~ 
20 days. In the context of a density jump, it is interesting to note that the inferred pre-jump 
circumburst density is se veral orders of magnitude lower than the values typically inferred for 
long-duration GRBs (see ISoderberg et al. Il2006d for a recent compilation). Moreover, radio 
observations of local (optically-selected) core-collapse SNe show similar flux modulations 



(factor of 2 to 3) attribu t ed to abrupt variations in the circumstellar medium (IWeiler et al. 
199ll : iRvder et aZ.ll2004j ; ISoderberg et al. Il2006bh . 



Next, the temporal and spectral evolution of the optical afterglow suggests the contribu- 
tion from a supernova component with peak luminosity and light-curve shape comparable to 
SN 1998bw. Given the temporal coincidence of the radio flare and X-ray flattening with the 
SN peak, XRF 050416a highlights the need for full spectral coverage in late-time GRB/XRF- 
SN searches. This is illustrated by the optical rebrig htening observed fo r XRF 030723 at ~ 15 
days and interpretted as a thermal SN component (IFynbo et al. 1120041 ) while the X-ray and 
radio dat a show similar rebrig htenings at late-time suggesting a CSM density jump or energy 
injection (IButler et al. II2005I . Soderberg et al. in prep). 



With the addition of XRF 050416a, there are three XRFs (020903, 060218 and now 
050416a) with spectroscopic redshifts observed in association with supernovae with peak 
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luminosities varying ; by up to a factor of 10 compared to SN 1998bw (ISoderberg et al. 112005b 



Bersier et al. 1120061 ; iPian et al. II2006I ). However, deep HST observations of XRF 040701 at 



0.21 revealed that an y associated SN was at le ast a factor of 10 (and likely > 100) times 



fainter than SN1998bw (jSoderberg et al. I l2005bl ). It is therefore clear that most XRFs 



produce Nickel-rich supernova explosions, but that there is a significant dispersion in the 
peak luminosities of XRF-assoc i ated SNe. A simi l ar resu lt is found for GRB-associated SNe 
( jZeh. Klose fc Hartmannl |2004| ; ISoderberg et al. I l2006dl ) and further strengthens the idea 
that GRBs and XRFs are intimately related. 

Finally, we address XRF 050416a within the framework of off-axis and dirty fireball 
models for XRFs. As discussed in §3.31 the evolution of the early afterglow indicates 
that the ejecta are being viewed along the collimation a xis. Similarly, on-axis v iewing an- 
gles are inferred fro m aftergl ow studies of XRFs 020903 (ISoderberg et al. Il2004af). 0502 15B 

Levan et al. Il2006h . 050406 faomano et al. Il2006h . and 060218 rtCanipana et al. koodlSoderberg et al. 
2006a ). Unification models in which X-r ay flashes are understood as typical GRBs viewed 
away from the burst collimation axis (e.g. Granot et al. 12002 ; Yamazaki. Yonetoku fc Nakamura 
20031 ) are therefore inconsistent with the observations for these five events. T o date, such 



models are only consistent with the observations of one X-ray flas h: XRF 030723 (IFynbo et al. 
2004! ; iButler et a/Tlbood : iGranot. Ramirez-Ruiz fe Perna lliool ). 



While the lack of high energy photons prevents a direct constraint on the Lorentz factor 
of the burst, our detailed afterglow modeling allows us to place a lower limit on the initial 
bulk Lorentz factor, T > 110, by extrapolating V oc t~ 3 ^ 8 back to the first XRT observation 
at ~ 10 sec. This value is somewh at lower than the typical values inferred for cosmological 
GRBs (IPanaitescu fc Kumarl 120021 ) and consistent with the upper range of values predicted 
for a dirty fireball. Give n that some dirty fireball models pre dict XRFs to have wider jets 
than typical GRBs (e.g.. IZhang. Woosley fc MacFadyenl 120031 ) . it is interesting to note that 
there are no X-ray flashes for which an achromatic jet break has been observed. This may 
suggest that their ejecta are not strongly collimated. In fact, detailed radio observations of 
XRFs 0209 03 and 060218 spanning > 100 days imply a quasi-spherical ejecta geometry in 
both cases (ISoderberg et al. Il2004al 12006a ). 



While the basic properties of X-ray flashes (redshifts, hosts, isotropic 7-ray energies) are 
now available for several events, it is clear that broadband afterglow observations are required 
for a complete understanding of the burst properties. With the addition of XRF 050416a 
to the existing sample of only a few well-studied events, we continue to work towards a 
systematic comparison of the global characteristics of XRFs and GRBs. While limited, the 
current sample suggests that the two classes share similar ejecta energies, associated SNe, 
and viewing angles. Further studies of XRF afterglows will be used to confirm whether the 
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two classes differ in their collimation angles. 
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Table 1. Ground-based Optical and NIR Observations of XRF 050416a 



Date Obs 




Telescope 


Filter magnitude 13 


(UT) 


(days) 







2005 


April 


16 


4641 





0025 


Palomar 60-inch 


/ 


18 


82 


± 





11 


2005 


April 


16 


4659 





0043 


Palomar 60-inch 


/ 


18 


86 


± 





11 


2005 


April 


16 


4677 





0061 


Palomar 60-inch 


/ 


19 


16 


± 





13 


2005 


April 


16 


4696 





0080 


Palomar 60-inch 


I 


19 


35 


± 





23 


2005 


Apri! 


16 


4692 





0076 


Palomar 200-inch 


K s 


16 


37 


± 





21 


2005 


Apri! 


16 


4714 





0098 


Palomar 60-inch 


I 


19 


01 


± 





12 


2005 


Apri 


16 


4731 





0116 


Palomar 60-inch 


I 


19 


53 


± 





18 


2005 


Apri! 


16 


4772 





0156 


Palomar 60-inch 


z' 


19 


16 


± 





27 c 


2005 


Apri 


16 


4763 





0147 


Palomar 200-inch 


K s 


16 


64 


± 





20 


2005 


April 


16 


4833 





0217 


Palomar 60-inch 


I 


19 


09 


± 





20 


2005 


April 


16 


4828 





0212 


Palomar 200-inch 


K s 


16 


79 


± 





20 


2005 


April 


16 


4880 





0264 


Palomar 60-inch 


I 


19 


34 


± 





23 


2005 


Apri 


16 


4895 





0279 


Palomar 200-inch 


K s 


16 


89 


± 





24 


2005 


Apri 


16 


4960 





0344 


Palomar 200-inch 


K s 


16 


99 


± 





25 


2005 


Apri! 


16 


5026 





0410 


Palomar 200-inch 


K s 


16 


98 


± 





29 


2005 


Apri 


16 


5089 





0473 


ANU 


2.3-meter 


R 


20 


92 


± 





05 


2005 


Apri! 


16 


5502 





0886 


ANU 


2.3-meter 


R 


21 


18 


± 





12 


2005 


April 


16 


5543 





0927 


ANU 


2.3-meter 


R 


21 


40 


± 





15 


2005 


April 


16 


5585 





0969 


ANU 


2.3-meter 


R 


21 


16 


± 





08 


2005 


April 


16 


5626 





1010 


ANU 


2.3-meter 


R 


21 


43 


± 





11 


2005 


April 


16 


5669 





1053 


ANU 


2.3-meter 


R 


21 


44 


± 





26 


2005 


Apri 


16 


5702 





1086 


ANU 


2.3-meter 


R 


21 


39 


± 





16 


2005 


Apri 


16 


5744 





1128 


ANU 


2.3-meter 


R 


21 


37 


± 





16 


2005 


Apri! 


16 


5779 





1163 


ANU 


2.3-meter 


R 


21 


18 


± 





17 


2005 


Apri 


16 


5816 





1200 


ANU 


2.3-meter 


R 


21 


35 


± 





16 


2005 


Apri! 


16 


5850 





1234 


ANU 


2.3-meter 


R 


21 


36 


± 





10 


2005 


April 


16 


5893 





1277 


ANU 


2.3-meter 


R 


21 


44 


± 





27 


2005 


April 


16 


5924 





1308 


ANU 


2.3-meter 


R 


21 


46 


± 





22 


2005 


April 


16 


5975 





1359 


ANU 


2.3-meter 


R 


21 


63 


± 





34 


2005 


April 


16 


6092 





1476 


ANU 


2.3-meter 


R 


21 


71 


± 





30 
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Table 1 — Continued 



Date Obs 




Telescope 


Filter 


magnitude 13 


(UT) 


(days) 








2005 April 16.6187 


0.1571 


ANU 2.3-meter 


R 


21.87 ±0.46 


2005 April 17.3895 


0.9279 


Palomar 60-inch 


R 


22.62 ±0.30 


2005 April 17.4006 


0.9390 


Palomar 60-inch 


z' 


< 20.38 c 


2005 April 17.4087 


0.9467 


Palomar 60-inch 


I 


22.05 ±0.30 


2005 April 18.2884 


1.8264 


Palomar 60-inch 


R 


< 23.46 


2005 April 18.4031 


1.9411 


Palomar 60-inch 


I 


< 22.83 


2005 April 23.3511 


6.8896 


Palomar 200-inch 


K s 


< 19.70 



a Days since explosion have been calculated for the mid-point of each 
exposure. 

b Magnitudes have not been corrected for extinction. Limits are given 
as 3a. 

C AB system. 
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Table 2. HST/ACS Observations of XRF 050416a 



Date Obs 


At 


Exp. Time 


Filter 


HST mag a 


Johnson mag b 


(UT) 


(days) 


(sec) 




(AB) 


(Vega) 


2005 May 23.38 


36.92 


3282 


F775W 


24.35 ±0.02 


23.82 ±0.02 


2005 May 23.46 


37.00 


3430 


F850LP 


23.83 ±0.03 




2005 June 16.17 


60.71 


3986 


F775W 


25.88 ±0.10 


25.36 ±0.10 


2005 July 11.15 


85.69 


4224 


F775W 


26.44 ±0.22 


25.92 ±0.22 


2005 November 21.06 


218.60 


4224 


F850LP 






2005 November 21.19 


218.73 


4224 


F775W 







a AB system magnitudes in the HST filters given in column 4. Photometry was done 
on residual images (see §2.2[) . We have assumed the source flux to be negligible in the 
final (template) epoch. Magnitudes have not been corrected for extinction. 

b We convert the F775W magnitudes in column 5 to Johnson /-band (Vega system) as 
described in §2.21 Magnitudes have not been corrected for extinction. 
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Table 3. Spectroscopic lines for XRF 050416a 



Line 


A (rest) 


A (observed) a 


Redshift 


Flux b 




(A) 


(A) 




(xl0~ 17 erg cm" 2 s _1 ) 


[Oil] 


3728.38 


6162.08 


0.6528 


9.6 ±0.4 


H 7 


4341.72 


7174.64 


0.6525 


1.2 ±0.2 


H/3 


4862.72 


8037.14 


0.6528 


3.7 ±0.5 


[OIII] 


4960.30 


8198.79 


0.6529 


2.5 ±0.3 


[OIII] 


5008.24 


8277.57 


0.6528 


8.4 ±0.4 



a Observed wavelengths have been corrected to vacuum. 
b Flux values have not been corrected for Galactic extinction 
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Table 4. Radio Observations of XRF 050416a 



Date Obs 


At 


^i/,1.43 GHz 


-Fi/,4.86 GHz 


F U ,8AG GHz 


(UT) 


(days) 


(^Jy) 


(^Jy) 


(^Jy) 


2005 April 16.49 


0.026 






20 ±51 


2005 April 22.03 


5.57 




200 ± 46 


101 ±34 


2005 April 28.29 


11.83 




188 ±42 


132 ±30 


2005 May 1.35 


14.89 




201 ±43 




2005 May 31.10 


44.64 






431 ±46 


2005 June 3.99 


48.53 




585 ± 48 




2005 June 8.03 


52.57 


420 ± 190 


562 ±51 


398 ± 33 


2005 June 16.97 
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189 ±36 
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41 ±47 
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156.44 




0±54 
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33 ±47 
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a All errors are given as la (rms). 




Fig. 1. — Discovery image of the afterglow of XRF 050416a. We began observing the field of 
XRF 050416a with the Palomar 60-inch telescope in the /—band about 2.5 minutes after the 
burst. We discovered a fading source within the 3.0 arcmin (radius) Swift/BAT error circle. 
Subsequent localizations of the XRT (circle) and UVOT afterglow positions were shown to 
be coincident with the P60 source. Our late-time HST images reveal a host galaxy coincident 
with the optical afterglow position. 
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Fig. 2.— Optical {BVRIz') and NIR (K s ) light-curves of the afterglow of XRF 050416a. 
We supplement our measurements from Table d] (b l ack symbols) with data from the GCNs 



(iLi et al. 11200 



2005 



| and 



)5; 



Price et al. 



Holland et al. 



2005; iQiu et al. 1120051 ; iToriil 120051 ; lYanagisawa. Toda fc Kawai 



( 20061 ) (grey symbols). Detections are shown as circles and upper 



limits as inverted triangles. R— band data from the ANU 2.3-meter have been binned for 
clarity. All data have been corrected for host galaxy extinction (see §2.11) . Between t ~ 0.01 
and 1 day the NIR/optical afterglow shows an average decay index of ttNiR/opt ~ —0.75. 
As described in §3.21 we find a reasonable fit to the spectral and temporal evolution of the 
broadband data with a standard afterglow model. Our HST measurements are shown as 
encircled dots; the implied steep spectral index suggests the contribution from an supernova 
component at t ~ 40 days ( §4.4j) . We sum the flux from a SN 1998bw-like supernova at z — 
0.6528 with our afterglow model to produce the final fits (grey solid lines). For comparison, 
we also show the afterglow model alone (dashed grey line). We note that this model does 
not apply to the very early afterglow evolution at t < 0.01 days. The dotted grey lines at 
early time represent the evolution of the NIR/optical assuming these bands track the X-ray 
evolution on this timescale. 
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Fig. 3. — The observed photometric spectrum for XRF 050416a is shown at t ~ 0.014 
(top) and 0.8 days (bottom). The K a — , z'—, I— and R— band data (black symbols) are 
from Table U while t he V, B, U, UVW1, UVM2 and UVW2 data (grey symbols) are from 
Holland et al. I (120061 ). Detections are shown as circles and upper limits as inverted triangles. 
The observed spectrum is fit by /^NiR/opt ~ —1-3 and —1.5 in the top and bottom panels, 
respectively (dotted line). We fit the dat a with an extinction model (solid line) which 
includes E(B — V) = 0.03 from the Galaxy (ISchlegel. Finkbeiner fc David Il998l ) and a rest- 
frame host galaxy extinction of E(B — V) rcst = 0.28 assuming an intrinsic spectral index of 
/^Nm/opt ~ —0.55 (dashed line). 



Q7 




Fig. A.— HST images of XRF 050416a were obtained with A CS in filters F775W and F850LP 
spanning t 37 to 219 days after the explosion. The host galaxy emission dominates that of 
the afterglow throughout this timescale. Relative astrometry between the P60 and the HST 
images provides an optical afterglow position accurate to 0.70 arcsec (circle, 2a). Residual 
images of the afterglow were produced through image subtraction techniques ( §2.21) . The 
position of the afterglow is offset by 0.02 ±0.02 arcsec with respect to the host galaxy center. 
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Fig. 5. — Optical spectrum of the host galaxy of XRF 050416a as observed with the Low 
Resolution Imaging Spectrograph on Keck I on 2005 June 6.3 UT. Several bright emission 
lines are detected, indicating that the host is a star-forming galaxy at redshift z = 0.6528 ± 
0.0002. As discussed in $5] we find that the host is actively forming stars with a rate of 
> 2 M Q yr _1 and an inferred metallicity of Z ~ 0.2 — 0.8 Z & . 
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Fig. 6. — Radio observations of XRF 050416a from the Very Large Array as listed in Table HI 
We have summed the late-time limits in each frequency for clarity. The radio evolution is 
characterized by an unusual flare at t ~ 40 days which rapidly fades below our detection 
limits. This radio flare occurs roughly on the same timescale as the observed X-ray flattening. 
(Figure [7]). We attribute both the radio flare and the X-ray flattening to an episode of late- 
time energy injection (refreshing shell or off-axis ejecta) or a large circumburst density jump. 
Our early afterglow model fit is shown (solid grey lines) and extrapolated to the epoch of the 
radio flare (dashed grey lines). The spectral index between 8.46 and 4.86 GHz is optically 
thin throughout our radio monitoring, including the late-time flare (lower right panel). 
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X-r a y ob servations of XRF 050416a from the Swift/XKT as reported by 



Mangano et al. I ( 120061 ). The X-ray light-curve is characterized by four phases: (1) a ini- 
tial steep decay, (2) a flattening between t ~ 0.005 — 0.01 days, (3) a subsequent decay, 
and (4) a second flattening between t 20 — 40 days. The final measurement at 70 days 
implies a rapid steepening following the second flattening. We attribute the late-time flat- 
tening to energy injection (refreshing shell or off-axis ejecta) or a large circumburst density 
jump. Overplotted is our afterglow model fit between 0.014 and 20 days (solid grey line) 
and extrapolated to late time (dashed grey line). We do not attempt to fit the data prior to 
0.01 days (dotted line) due to insufficient broadband data. 



